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Abstract 

Dynamical evolution of the quantum ground state (vacuum) is analyzed for time 
variant harmonic oscillators characterized by asymptotically constant frequency. The 
oscillatory density matrix in the asymptotic future is uniquely determined by a 
constant number of produced particles, independent of other details of transient 
behavior at intermediate times. Time average over one oscillation period yields a 
classical, in some cases even an almost thermal behavior. In an analytically soluble 
model the created particle number obeys the Planck distribution in a parameter 
limit. This suggests a new way of understanding the Gibbons-Hawking temperature 
in the de Sitter spacetime. 
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There are a number of situations in which a precise knowledge of evolved quan- 
tum state in the simple system of time variant harmonic oscillators is expected 
to deepen our understanding of fundamental problems in physics. A class of such 
problems of great importance is emergence of classical, in some cases even thermal 
behavior, resulting from time evolution of pure quantum ground state in a nontrivial 
spacetime. The central issue here is the behavior of quantum field in the curved 
spacetime, the quantum field being an infinite collection of harmonic oscillators with 
different frequencies, interacting each other in some prescribed way. One may men- 
tion as examples of such problems the thermal Hawking radiation from formation 
of black holes Jj]], and quantum dissipation and fluctuation in the inflationary de 
Sitter universe @. Although differences in indivisual cases are important, there are 
some common features that deserve special attention. In the classical approximation 
to gravity the spacetime in the examples listed is highly curved or is endowed with 
event horizon. In an idealization that extracts essential features of the underlying 
physics, one replaces the background spacetime by a simpler time evolving space. 
For instance, the black hole formation may be modeled by a space bounded by mov- 
ing mirror [[J that acts as the source of red-shift just as in the gravitational potential 
of the black hole. Investigation of the behavior of quantum fields in such a moving 
mirror spacetime has led to Hawking radiation by taking the mirror trajectory cor- 
responding to the black hole formation. This picture is ragarded as an essence of 
the thermal behavior of the black hole spacetime. 

In this paper we would like to lay down foundation of these and other related 
problems by solving the most basic and the simplest problem of this nature: time 
evolution of the quantum ground state for time variant harmonic oscillator. Time 
variation of harmonic oscillator frequency is taken arbitrary here, provided that the 
initial and the final frequencies approach constant values. The same problem has 
undoubtedly been examined many times in many areas of physics, but to our knowl- 
edge the final quantum wave function or its density matrix has never been worked 
out with sufficient generality. For instance, Parker and Zeldovich M only discussed 
rate of particle production in an expanding universe. Our construction of the state 
vector clarifies more details; we shall be able to write down the quantum density 
matrix itself in terms of a few parameters in such a way that detailed knowledge of 
time evolution at intermediate times is not needed. Only with such an apparatus 
one may discuss the delicate issue of transition from the quantum to the classical 
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behavior with some kind of coarse graining. We consider in this paper the most 
straightforward coarse graining: time average over one oscillation period. This av- 
erage is also related to a coarse graining in the phase space variable of the Wigner 
function. Even this simple time average yields a new insight into the fundamental 
problem set forth; how the thermal behavior emerges from a quantum state. Our 
emphasis is on dynamical time evolution, rather than on geometry such as the event 
horizon discussed in many of the recent literature 0,0], 0. 

To illustrate our method of analysis we also briefly sketch an exactly soluble 
model of time variation. When suitably interpretated in a parameter limit, this 
model describes evolution of the quantum ground state in the de Sitter spacetime 
relevant to the inflationary universe. The final state in this limit is a nearly thermal 
one, nearly because the distribution of produced particles is exactly Planckian, but 
the time averaged density matrix has a deviation from the thermal one. Our model 
suggests a new way of understanding the Gibbons-Hawking temperature |J, but 
further work is needed to conclude more definitely. 

Quite apart from our motivation described above, our explicit solution to the 
problem of time variant harmonic oscillator is expected to be useful in solving other 
fundamental and practical problems in many areas of physics. With this in mind we 
stress the technical aspect of our approach in the main body of the paper. In the end 
of this paper we come back to physical applications of the mathematical analysis. 

Consider the generalized harmonic oscillator given by the following Hamiltonian: 

H=\p* + \u J \t) q \ (1) 

The variant oscillator frequency u(t) may have an arbitrary time dependence pro- 
vided that it has a definite limit in both the asymptotic past and the future: 
u(t) — > ujq as t — > — oo , u(t) — ► u) as t — > oo . For a general u(t) the wave 
function at any finite time is difficult to write down in a useful form. But as we shall 
explicitly construct in the following, the asymptotic form of the transition amplitude 
can be worked out with only a few parameters by using the path integral formula. 
We shall slightly extend the original Gelfand-Yaglom method |7| by simultaneously 
solving two boundary value problems, with the Cauchy data given either at the given 
initial time tj or at the given final time tp. 

The starting point of this construction is the expression for the transition ampli- 
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tude when the discretized time slicing with a step e = (t F — tj)/(N + 1) is employed: 

1 A dq t 



f 1 N 

q F ; t F | g 7 ; tj ) = lim / -== JJ 



exp[ | £ fti^fc - V g + ± (i^g 2 + tf"g 2 ) ] . (2) 

The matrix X y has nonvanishing elements only along the diagonal and the next 
off-diagonal lines: K" = 2 — t 2 uof , K l,l+1 = K l+1 < l = — 1 . Furthermore K FF = 
2 — e 2 cj^ +1 , K 11 = 2 — e 2 a; 2 , c = ( q z , , • • • , , q F ) . Integration with respect 
to the oscillator coordinates g, yields the transition amplitude in terms of the inverse 
matrix K~ l : Ki} = N ~' , for i > j . Df' is the determinant made out of a 
part of the matrix K of dimension i, in the upper left (u) or in the lower right (d) 
corner. They obey recursion relation with smaller dimensions: = K %% D^\ — Di_ 2 
for both £)" ' d . 

In the continuum limit of e — > the recursion relation combined with the end 
point value of D u > d leads to the Gelfand-Yaglom equation with a definite boundary 
condition. Thus with Di = eDj 

(^ + u 2 (t))D(t) = 0. (3) 



The boundary conditions are 

D"{t I ) = 0, ±D"(t z ) = l t (4) 
D d (t F ) = 0, ±D d (t F ) = -l. (5) 

The continuum limit of the transition amplitude is then given by 

/ c \ 1/2 i 

( q F ; t F | q r ; tj ) = [—) exp[ - (Aq] + Bq 2 F - 2Cqiq F ) ] , 

A = -4- ^D d (tj), B = -4- \nD u (t F ), C=^— = ^ — .(6) 
dt T V /; dt F K Fl ' D d {ti) D u {t F ) y 1 

In the rest of this paper we shall be concerned with evolution of the ground state. 
The evolved ground state at time t F is obtained, either by multiplying to Eq.^ the 
initial ground state wave function and integrating with respest to qi, or by directly 
taking the analytically continued limit of itj — > — oo in Eq.|6]. 

One may consider an analogue potential model of zero energy under the potential 
— u 2 (x) by replacing the time t by a space coordinate x. In the far left (x < 0) 
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region the potential approaches — u 2 , while in the far right — uj 2 . The Gelfand- 
Yaglom boundary condition is very specific, for instance, D u (x) must have the form, 
sin(a;o(x — ti))/u> in the far left near x = tj. 

In solving the Gelfand-Yaglom problem it would be of great value if one can de- 
couple the detailed transient behavior at any finite time and concentrates on asymp- 
totic past and future regions. This approach is useful unless one asks minute details 
at intermediate times. We thus parametrize the asymptotic form of solution respect- 
ing the specified boundary condition. Let us discuss D u (t) in some detail. It is real 
by definition and is given in terms of two linearly independent solutions, ip and ip* 
to the differential equation Eq.|3|: 

D u (t) = if; + i)* -> — — ( e^*-*^ - e" ) , at t ~ h ~ - oo , (7) 

i)(t) -> ai e iujt+iSl +a 2 e- iujt - i&2 , at t ~ t F ~ oo , (8) 

where ai is defined as both real and positive, and Sj are real phases. Of course, with 
given frequency variation u(t) these parameters Oj , Si are determined by solving the 
differential equation. Without resort to explicit form of the solution we shall be able 
to derive certain useful results by considering general properties of these quantities. 
We immediately note that a 2 — a 2 = , verified from constancy of the Wronskian, 
and the phases Si depend on the initial time like S\ ~ — u tj , S 2 ~ uJoti . 

Similarly one derives the asymptotic form of D d . Starting from 
D d (t) ~ e laJ (* F- *) + (c.c.) near tp, one defines at t ~ tj ~ — oo 

D d (t) — > a ie - iuJot ' 1 ' 51 + a 2 e itJ0t+fS2 + (c.c.) . (9) 

In terms of these parameters one finds that the quantities in the transition amplitude 
(Qf )tp \ qi ;ti) in Eq.| is given by 

A = Uq COt(u)otl + (fl) , B = UJ COt(u)tF + <Pf) , 

C=[2\ a x e lSl + a 2 e i&2 \ sm(ujt F + <p F ) }~ 1 = [ 2\ a x / Sl + a 2 e f52 \ sin(w t/ + <Pi) ]" 1 , (10) 
where (p F is independent of t F and 

ai sin 5i + a 2 sin S 2 

-cot(p F = — — . 11 

ai cosoi + a 2 coso 2 

A similar relation holds for (ft with and Si replaced by respective tilde quantities. 
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One may take the tj — ► zoo limit in order to extract the wave function of the 
evolved vacuum state. Dependence on the initial time tj in Eq.|I], besides the explicit 
one, is through Si as already mentioned. It is convenient to define a quantity A by 

d = — —= cosh A, a 2 = — -= sinhA. (12) 

In the analogue potential model tanh 2 (A) = (a 2 /ai) 2 is the fraction of left-mover 
relative to right-mover in the far right region. A straightforward computation shows 
that the wave function of the evolved vacuum state is given by 

( q F t F | ) = ( irrelevant phase factor ) x 

/ 7T . , / . x , / « % n \ _1 ^ 4 r ^9 1 + i sinh(2A) sin F , , . 

- cosh 2A + sinh 2A cosF exp - - ql — —-- . . , 13 

V u V V ; K J 1 ) [ 2 HF cosh(2A) + sinh(2A) cosF J ' v ; 

where the oscillaroty factor is given by 

F = 2ut F + ( phase of a ie i5l+iuJotl ) + ( phase of a 2 e iS ' 2 - iuJotl ) . (14) 

The Wigner transform of the density matrix, namely the Wigner function, takes a 
particularly simple form. It is calculated as 

fw(q,P;t F ) = 2exp[-- cosh(2A) (p 2 + u 2 q 2 ) - - sinh(2A) (p 2 + u 2 q 2 ) cos(F + X ) ] , 

OJ OJ 

(15) 

where tan x = 2ojqp/ (oj 2 q 2 — p 2 ) . 

A remarkable feature of the general formula just derived is that the wave function 
and the Wigner function depend on the detailed transient behavior at intermediate 
times only via the quantity A. The other phase factor is contained in F Eq.[T4] 
together with the final time tp. Furthermore the role of the Hamiltonian H is 
distinct in EqjTjJ all points on the energy surface of fixed H in the phase space 
(q , p) run periodically in the same range of the Wigner function, the exponent varying 
sinusoidally within — 2e 2A H/oj ~ — 2e~ 2A H/oj . Compare this with the value for the 
quantum ground state, — 2H/oj, and that for the thermal state, — — tanh(/3tt;/2) . 
Time variation of the potential raises the possibility of particle production in the 
evolved state from the initial vacuum, which is indeed what takes place in this case. 
The quantity A is related to the final number of produced particles as follows. By 
using the standard formula of the creation and the annihilation operators a w = 
— ? Vf 1 + ~j2^Pi one com P u tes, with the wave function Eq.|13"l, ( a^a^ ) (t F ) = 
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sinh 2 A, which is also equal to Aujqujci 2 .. The produced number of particles is thus 
independent of time tp asymptotically: no oscillation occurs for this measurable 
quantity. 

Rapid oscillation of the Wigner function in the time scale of 1/u naturally raises 

the question as to its time averaged behavior. Certainly one loses important quantum 

information by the procedure of time average, but some kind of coarse graining is 

of prime importance in statistical interpretation. The time average of the derived 

Wigner function Eq.[T5| over one oscillation period is 

1 f 47r , ^ „ , . r cosh2A . 9 9 9n . T r sinh2A . 9 9 9 . , . , 

— / dFf w (q,p;t F ) = 2ex V [ (p 2 + uj V) Ro (p 2 + u 2 q 2 )}, (16) 

An Jo uj uj 

where Iq(x) is the modified Bessel function. This time-averaged Wigner function is 

diagonal in the Hamiltonian basis, and in the large argument limit it reads as 

/2tt \~ 1/2 1 

2 — sinh(2A) (p 2 + u 2 q 2 ) exp[ - - exp(2A) {p 2 + uj 2 q 2 ) ] , (17) 
\ uj J UJ 

for p 2 + uj 2 q 2 ^> 1. The density matrix p corresponding to Eq.pl] does not obey the 

relation, p 2 = p, required for a quantum state, due to the average procedure. The 

time average over the period can be shown to be equivalent to a coarse graining over 

the phase space region of variables ( q , p) of the Wigner function. Namely, from the 

X dependence of Eq.|T^ one finds the identical coarse grained Wigner function EqirBj 

by taking the average in the range of \Ap\ = uj\Aq\ = \JlH on the fixed energy (H) 

ellipse. This exactly corresponds to the classical phase space motion over one period 

for the fixed frequency uj. 

In physical applications of the mathematical analysis above it often becomes 
necessary to extend to a system of many oscillators, infinitely many in the case 
of quantum field theory, with further complication of interaction among them in 
any realistic problems. Even if this coupled system follows a universal change of 
oscillation frequency given by a single function of time, one must be careful in taking 
the time average. The simple average formula given above is valid for time integration 
over time interval > 2it/uj. Thus the low frequency part must be treated separately 
if the corresponding period is longer than the time interval considered. The infrared 
component of the quantum field thus needs a special care. 

To gain more insight it is instructive to explicitly work out the following model 
of time variation: 

uj 2 (t)=uj 2 i Kt , a>0, k>0. (18) 
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The oscillation frequency monotonically changes with a time scale 1/k like uj(t) = 
Uq — > ^/acuo = to as time evolves. It turns out that this model is analytically soluble 
at any finite time ||. The solution to the Gelfand-Yaglom problem is expressed by 
the hypergeometric function. For instance, 

D u (t) = i^e*"* F( i^(VE + 1) , - z-(VS - 1) , 1 + i— ; - e Kt ) + (c.c.) , (19) 

ZILUq k k k 

with F( a , (3 , 7 ; z ) the hypergeometic function. For £/ — > — 00 , e = O (e Ktl ) , hence 
e may be ignored. Using the analytic continuation formula of the hypergeometric 
function, one finds for the parameters in Eq.[8] that 

a ie * = ^-e~ *** 2^(^ + D B( ^ , ^ ,/a ) B( i + ,J + <^ Va ) , 

47TK K K 2 « 2 K 

a 2 e" * = J- e- 2" ^<^> S( , - ^ ^ ) S( - + , - - ) , (20) 

A7XK K K 2 K 2 ft 

with S(a , /?) the Euler's beta function. The crucial quantity A = arctanh (02/01) is 
then given by 

tanhA= |sin ^ (7r / ( "-^ )l . (21) 
smh( Ti (u + ojq)I k ) 

From this one obtains for the created particle number 



, . , 1 cosh 2z(uj — u Q ) — 1 

< N » ' = 2 sinh^sinh^ • (22) 



Let us take a limit of ujq — > 00 such that the final u; = y/auo is kept finite. Shortly 
afterwards we shall argue that this limit is relevant to the de Sitter expansion. The 
limit gives tanh A — > e~ 2 ^ U> / K ; and the produced number of particles ( ) — > 
e 4nJ/K_ 1 ■ This is exactly of Planckian form with the temperature: T = ^ . 

Let us now discuss an application of the analytic model just solved. In the a — > 
limit there exists an intermediate epoch in which e Kt ^> 1 but ae Kt <C 1, namely 
<C — In a/ft. In this epoch uj{t) oc e~ Kt l 2 . This behavior may be considered 
as an analogue model of de Sitter inflation. Suppose that all modes of quantum 
fields universally follow this exponential decrease of frequency, as occurs for massless 
fields. The parameter k is then interpretated as the Hubble rate of expansion: 
H = k/2. (Strictly, modes of large physical wavelength > l/H do not experience 
the exponential frequency decrease, and must be treated separately.) After the epoch 
of de Sitter phase the space expansion ends and the universe settles to a mush slower 
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phase of expansion. If this final universe and also the initial one are approximated 
by the static universe, one may consider the final produced particle being created by 
the de Sitter expansion. Indeed, the resulting Planck distribution coincides with that 
of Gibbons-Hawking [0] , including the temperature relation: T = 3- . We find this 
interpretation worthy of further investigation, because in the original approach of 
Gibbons and Hawking the thermal radiation is observer dependent and its reality is 
not unambiguously established. On the other hand, in our model the time averaged 
density matrix is nearly, but not completely thermal. This and the infrared problem 
mentioned above must further be examined. 

In summary, we have given the general formula of the transition amplitude for 
time varying harmonic oscillator, and in particular the wave function of the evolved 
quantum ground state Eq.|13| and its Wigner function Eq.|15[ The resulting number of 
created particles approaches a constant value if both the initial and the final oscillator 
frequencies have well defined limits. We have also discussed how the thermal behavior 
may emerge from a quantum state via coarse graining. 
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